Active Particles Forced by an Asymmetric Dichotomous Angle Drive 
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We analyze the dynamics of particles in two dimensions with constant speed and a stochastic 
switching angle dynamics defined by a correlated dichotomous Markov process (telegraph noise) 
plus Gaussian white noise. We study various cases of the asymptotic diffusional motion of the 
particle which is characterized by the effective diffusion coefficient. Expressions for this coefficient 
are derived and discussed in dependence on the correlation time and the intensity of the noise. 
The situation with a given mean curvature is of special interest since a non-monotonic behavior of 
the effective diffusion coefficient as function of the noise intensity and correlation time is found. A 
timescale matching condition for maximal diffusion is formulated. 
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Introduction. We study self-propelled particles 
moving in two dimensions at a constant speed Vq. 
The time-dependent position r(t) of a particle follows 
from integrating its velocity v(t) — (v x (t),v y (t)) — 
vo (cos sin 0(f)) with time-dependent orientation 
4>(t). The orientation of the velocity vector at time t is 
governed by stochastic dynamics. We assume that <p(t) 
changes due to a constant torque superimposed by an 
unbiased dichotomous Markov process (DMP) £(i) which 
increases or decreases the local curvature of the particle's 
trajectory. In addition, a Gaussian white noise is present, 
corresponding to the thermal or environmental noise of 
the system. 

Physically, this dynamics is motivated as an ap- 
proximation to recently measured bimodal distributions 
P(A, r) of turning angles A = <fr(t + r) — <p(t) during 
time r, as observed in experiments with the zooplankton 
Daphnia [H-01, which is also able to sustain a constant 
mean speed over large time scales. By using a DMP, we 
approximate the bimodal structure by two delta peaks 
at (Q + A_)t and (f2 + A + )t, where and A + denote 
the DMP strokes and 51 is the additional torque-induced 
angular velocity, which we will simply term "torque" in 
what follows. The constant torque can be motivated by 
various biological realizations. On the one hand, there 
are typical swarming characteristics which can be intro- 
duced by an effective torque [l|, On the other hand, 
an asymmetric muscularity Q , an external magnetic field 
Q or corresponding asymmetric boundary conditions Q 
can lead to an effective torque as well. 

Thus, our system is fully described by the constant 
speed \v\ = vq, and the angle dynamics 



0(t) = n + C(t) + £(<). 



(i) 



is the Gaussian white noise with zero mean and noise 
intensity D^. 

The dichotomous Markov process £(i) is time- 
homogeneous and switches between the two values A + 
and A- with transition rates r + and r_ [8|,[9[. r + denotes 
the rate for changing from A + to A_ and r_ denotes the 
rate of passage from A_ to A + . The mean value of this 
process, i.e. (C(^)) = ( r -A+ + r +J 4_)/(r + +r_), is fixed 



20 



10 



-10 




-20 



20 



FIG. 1. (Color online) Spatial trajectories for dichotomous 
angular dynamics (see Eq.Q; details are shown in the inset) 
with = 0.01, vo = 1, vanishing torque O = 0, and the 
total duration ti = 10. Two parameters of the DMP are fixed 
to be A- = —2, r+ = 2, A+ varies, and r_ is defined by the 
vanishing mean of the DMP drive. 



at zero in what follows, since the corresponding mean 
value can always be incorporated into the term Q. 

At first, we derive an explicit expression for the effec- 
tive diffusion coefficient in our system which is defined as 
the long-time limit 



lim 

t— ^oo 



m)-r ) 2 ) 

At 



(2) 



We then discuss its dependence on the DMP parameters 
and on the external noise for the case without torque as 
well as for the case with non- vanishing torque Q. We 
observe a torque-induced non-monotonic behavior of the 
effective diffusion coefficient, similar to the one recently 
discussed in the context of a system which is driven by 
an Ornstein-Uhlenbeck process (OUP) [lfl], as well as to 
the peaked diffusion of spiral waves driven by a correlated 
random forcing [ll|. Finally, we study the DMP limits 
leading to shot noise and to Gaussian white noise [f3 |. 
The latter reproduces the well known result which was 
previously derived in [Io|, [H, [3 ■ 
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Analytical considerations of a DMP-driven 
agent. In Fig[TJ we present the spatial trajectories as 
obtained from simulations of our dynamics [Eq.Q] for 
vanishing torque, small noise intensity and different 
parameters of the dichotomous noise. The trajectories 
clearly show the dominant circular structures, as well as 
the strong influence of the parameters of the DMP on 
the particle's displacement. The effective diffusion coef- 
ficient of particles moving at a constant speed vq, can be 
written with the Taylor-Kubo relation as 

-| poo 2 poo 

D cS = -J {v(t)v{t + r))dr = f J (cos(A(r)))dr 

= f Re QT J°° e lA P(A, r)dAdr) , (3) 

where P(A,r) denotes the probability for an angle in- 
crement A during time t. Using a continuous-time gen- 



eralization of the classical persistent random walk, as it 
was studied in [l5j], we derive the probability to finish a 
dichotomous Markov step at a certain time with a cer- 
tain angle. This calculation can be explicitly done Q and 
leads to a coupled system of integral equations, which can 
be solved in an algebraic way by considering the Fourier- 
Laplace transform 

/oo />oo 
e lkA / e- ST P{A, r)drdA (4) 
- oo JO 

of the corresponding probability density. With P(k,s), 
Eq.Q can be rewritten as 

v 2 

Deff=yRe(P(* = l,s = 0)) ) (5) 

and we are finally able to derive the effective diffusion 
coefficient of our dynamics 
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n ^ (gg + r + + r_)(£> e (£> g + r + + r_) - (A + + 0)(A_ + Q)) 

"" ^ (D € (£»£ + r+ + r_) - + fi)(A_ + f!)) 2 + + n)(D £ + r_) + (A. + fi)(I> 4 + r+)) 2 

( A+ ;+y_- r - + n) ((A + + Q)(D t + r_) + (A. + + r + )) 



(D C (I>4 + r+ + r_) - (A+ + ft)(A_ + fl)) 2 + ((A+ + Q.)(D e + r_) + (A_ + ft)(£> ? + r+)) 5 



(6) 



Without loss of generality we assume that A + > > A_ and \A + \ > \A_\ which also implies r + > r_. For describing 
the DMP dynamics in a more intuitive way, we use the following parameters fill ] 

A = A+-A_, r c = - , 0<p = --t^<1. (7) 

r+ + r_ A + 

Now, A measures the strength of the process and r c is the correlation time of the DMP. The parameter p controls 
the asymmetry of the driving, the symmetric case corresponds to p — 1, and it tends to for strongly asymmetric 
strokes. 

Rewriting Eq.Q in terms of these new parameters, we get 



Dt + ) Di + 
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Discussion of the effective diffusion coefficient. 

Considering a vanishing torque = 0, we can discuss the 
symmetric and asymmetric limit of our DMP drive. A 
strong asymmetry leads with Eq.© to 

so that we receive a divergent expression for a vanish- 
ing intensity of the thermal noise £(£). Hence, the noise 
leads to the maintenance of the diffusive character of our 
system. This limit coincides with the result of the Gaus- 
sian white noise-driven angle discussed in [l7i j . In the 



limit p — > 0, the DMP-term £(i) becomes negligible for 
the simple reason that the torque A- vanishes. There- 
fore, only the Gaussian force drives the angle which 
results in Eq.©. 

The symmetric, torqueless case p = 1 on the other 
hand, reproduces the result 

,,2 D f + — 

which has been previously derived in [2j. Analyzing this 
expression leads to a noise intensity which maximizes 
the effective diffusion coefficient, namely the value which 
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FIG. 2. (Color online) Spatial trajectories for different r c val- 
ues and a non- vanishing torque Q = 1; the other parameters 
are p — 0.4, = 0.01, A = 2, vq = 1, and total duration 
t; = 10; the inset illustrates the effective diffusion coefficient 
D c fi versus correlation time r c for different asymmetry pa- 
rameters p. 
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FIG. 3. Effective diffusion coefficient versus noise intensity 
(left) and torque Q (right) for different asymmetry parameters 
p within the theory (cf. Eq.©) for the constant torque f2 = 
l(left), D ( = 0.01 (right), r c = 1, v = 1, A = 2. 



obeys + l/r c = A/2. The two dissipative time scales 
of the noises sources, i.e. = 1/-D{ and r c , relate to the 
torque of the DMP which turns the particle. 

A non-vanishing torque SI ^ changes the basic char- 
acteristics of our dynamics and induces a systematic 
change of the curvature of the path. Figure [5] shows 
spatial trajectories for different values of t c and illus- 
trates the influence of fi. The fact that this influence on 
-D e ff leads to a non-monotonous parameter dependence is 
demonstrated in Figl^inset) and [3] We recognize that 
D e ff as a function of r c , and SI possesses a well-defined 
maximum (a similar p dependence is not shown here). 

This can be understood in view of the corresponding 
behavior of the agent considered [see Fig(5]. Small and 
large correlation times lead to a curled structure where 
the system stays in a certain DMP state either too long 
or too short in order to perform a considerable spatial 
displacement. In the case of r c = 10, the particle per- 
forms a persistent circular motion with small curvature 
(SI + A_) interrupted by small spins with large curva- 



ture (SI + A + ) and the additive noise causes diffusion 
by shifting the centers of the circles stochastically. For 
the case with r c = 0.1, fast DMP-switches induce an er- 
ratic motion. In contrast to the torqueless situation, the 
non- vanishing SI reduces the displacement of the particle. 
On average, the motion follows again randomized circu- 
lar lines determined by the non-vanishing torque and the 
fast DMP strokes, whose mean influence disappears for 
fast switchings. Calculating the limits of large and small 
r c values analytically, results in a non-zero value of D c g 
in both cases due to the additive Gaussian noise. This 
property is seen in Fig J5] where both asymptotics tend to 
finite values. 

The optimal r c in between the two limits, i.e. the 
state of maximal diffusion, corresponds to a maximally 
stretched trajectory for given values A± and only the A- 
stroke can decrease the curvature. Thereby, it can induce 
longer excursions which become maximal if the mean 
waiting time l/r_ matches the time which the angle <p 
needs to rotate over half of a circle during the A_ stroke. 
Hence, the parameters have to obey \A^ + Sl|/r_ = it 
and for the notation introduced in Eq. (0 follows 
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This result is in good agreement with the peaks in FigJ2] 
for small p values. It fails for larger p, where the influence 
of the A + stroke is not negligible. 

The behavior of D c s as a function of SI [Figj3|right)] 
shows a peak due to similar reasons. The case SI = |^4-| 
causes straight paths within the corresponding DMP 
mode and will therefore enhance the spread. The peak 
in the dependence of D e g on the noise intensity in 
FigEIleft) shifts for growing p to smaller noise values. 
Such a peak was already reported for similar dynamics 
but in the absence of a DMP [ID, [l4[ ■ It is in agreement 
with our previous discussion of Eq.® that the effect of 
the DMP strokes disappear if p —> . 

Taking the derivative of Eq. ([5]) with respect to r c leads, 
after some straightforward calculations, to a lengthy ana- 
lytical result. t c values which maximize the diffusion co- 
efficient are presented in Fig[4j It shows a perfect agree- 
ment with simulation results and with the peaks in FigJ5] 
The rough approximation given by Eq. (|1 1[) turns out to 
be rather good. The peaks in Figd] occur because of the 
mentioned rectilinear motion for ft = \A- \ = pA/(l +p), 
since an increase of r c also increases the duration of the 
rectilinear motion and enlarges consequently D c g. That 
is why the peak for p = 1 is shifted to infinitely large 
correlation times in FigJSJ inset). 

The white shot noise limit of the DMP drive can be 
found by considering the limits A+ — > oo and r + — > oo 
while the ratio A + /r + = —A-/r- = w holds constant 
The corresponding locomotion of the agent con- 
sists of a circular motion with mean curvature (Q + A_ ) 
interrupted by infinitely fast turnings of the angle, in- 
duced by the short, large strokes A + . The autocor- 
relation function of Q(t) in this shot noise limit reads 
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FIG. 4. (Color online) r c value which maximizes the effective 
diffusion coefficient versus asymmetry parameter p for differ- 
ent mean torques f2; the results of the full theory in Eq.® 
(solid lines) and for the approximative expression in Eq. (|ll[) 
(dashed lines) are shown for the noise intensity = 0.01 
and the stroke strength A — 2. 



(C(i)C(* + t)) = (-A_w/T c )exp(-\T\/T c ) and therefore 
implies the noise intensity = — A-w, while < 
holds. Rewriting of Eq.Q within the mentioned limits 
leads to the effective diffusion coefficient 
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which implies < \A- 

The white Gaussian limit of the DMP drive can be 
derived b y c onsidering the limit w — > while is hold 
constant [12[- Doing so in Eq. (jT2"j) . we find 
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If we introduce a total noise intensity D = + Dr, this 
expression coincides with the result in [l(], Il3l [l4J | for 
an agent under influence of Gaussian white noise with 
intensity D. Maximal diffusion at the value Z?^jj uss = 
Vq/(4H) is obtained for the total noise intensity D = 0. 
Thus, the resonance occurs where D, i.e. the angular 
correlation decay rate in the case of a Gaussian white 
angle drive, equals the effective torque, as it is likewise 
the case in Eq. (fT0|) with the additional decay rate of the 
correlation within the DMP drive. 

Conclusion. We have discussed exact results for the 
effective diffusion coefficient of a particle moving at a 
constant speed under influence of a constant torque, di- 
chotomous angular Markov noise and additional direc- 
tional Gaussian perturbations. The results help to un- 
derstand the behavior of our system in a qualitative and 
quantitative way. They clarify the role of the asymme- 
try and of an additional torque in the DMP-driven angle 
dynamics. The strongly peaked bimodal angular proba- 
bility distribution, which we have assumed in our model, 
is of course an enormous simplification of the ones found 
in real biological systems. But in view of the bulk of 
works discussing symmetric angular distributions 0, Il7r - 
[2?jj . it seems to be reasonable to discuss the influence of 
a certain asymmetry as well. Since the diffusion coef- 
ficient is one of the most easily accessible quantities in 
experiments, we hope that this work not only fills a gap 
in our general theoretical understanding of self-propelled 
agents, but will also stimulate corresponding experimen- 
tal studies. 
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